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GENERAL  VARIATIONAL  METHODS  FOR 
WAVES  IN  ELASTIC  COMPOSITES1 

by 

S.  Nemat-Nasser3 

ABSTRACT 

General  variational  theorems  in  which  the  displacement,  the 
stress,  and  the  strain  in  one  case,  and  the  displacement  and  the  stress  in 
another  case,  are  given  independent  variations,  and  which  include  appro¬ 
priate  general  boundary  and  discontinuity  conditions,  are  developed  with 
a  view  toward  the  application  to  harmonic  waves  in  elastic  composites  with 
periodic  structures.  The  one- dimensional  case  is  first  developed  in  detail, 
and  in  order  to  demonstrate  the  effectiveness  of  the  results,  especially 
their  accuracy  in  providing  the  dispersion  curve,  waves  propagating  normal 
to  layers  in  a  layered  composite  are  discussed,  and  numerical  results  are 
presented;  see  Tables  I  and  II.  Then  the  general  three-dimensional  case 
is  considered,  and  the  results  are  applied  to  waves  propagating  normal  to 
the  fibers  in  fiber- reinforced  composite. 


1  This  work  was  partly  completed  while  the  author  was  at  the  University  of 
California,  San  Diego,  La  Jolla,  California,  as  a  consultant  to  Grant 
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Professor,  The  Technological  Institute,  Northwestern  University, 
Evanston,  Illinois., 


I. 


INTRODUCTION 


Because  of  their  dispersive  and  other  desirable  effects,  composite 
materials  have  become  now  an  important  ingredient  in  many  aerospace 
structures..  The  propagation  of  elastic  waves  in  such  materials  has,  there¬ 
fore,  been  discussed  by  a  number  of  writers  using  different  theories.  In 
[1,  2],  for  example,  a  two- term  Taylor  expansion  of  the  displacement  field 
within  each  layer  in  a  layered  composite  is  considered  together  with  a 
certain  smoothing  process,  to  arrive  at  an  approximate  set  of  equations 
which  resemble  those  for  a  homogeneous  continuum  with  microstructure. 
Bedfoid  and  Stern  [3],  on  the  other  hand,  directly  consider  a  mixture 
theory  to  characterize  the  composite,  and  proceed  to  calculate  the  coupling 
coefficients  in  the  corresponding  constitutive  relations,  by  considering 
simple  static  elasticity  problems;  these  authors  confine  their  analysis  to  a 
special  case  where  the  interacting  body  forces  are  assumed  to  be  propor¬ 
tional  to  the  relative  displacement  of  the  constituents.  In  a  more  recent 
paper,  Hegemier  and  Navfeh  [4]  have  used  an  asymptotic  approach  and,  for 
waves  propagating  normal  to  the  layers  in  a  layered  composite,  have  derived 
systematically  mixture-type  field  equations  directly,  which,  for  the  har¬ 
monic  wave,  gives  the  exact  dispersion  relation,  and  which  can  be  used 
(as  these  authors  do)  to  study  the  trans ient  waves  in  such  composites. 

For  harmonic  waves  in  a  composite  with  a  periodic  structure,  one 
may  employ  a  variational  approach  This  can  particularly  become  a  very 
effective  tool,  if  one  uses  a  variational  statement  in  which  not  only  the 
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displacement,  but  also  the  stress  field  is  given  independent  variation. 
Moreover,  by  permitting  discontinuity  in  the  displacement  and  the  stress 
test  functions,  one  can  expect  a  more  accurate  reproduction  of  the  local 
variation  in  the  displacement  and  the  stress  fields  within  and  across  the 
constituent  materials.  An  example  of  such  a  calculation  can  be  found  in  a 
thesis  by  Wheeler  [5]  and  in  a  recent  article  by  Kohn,  et  al.  [6],  where 
the  classical  theorem  of  stationary  potential  energy  which  leads  to  the 
well-known  Rayleigh  quotient  for  the  eigenfrequency,  is  used.  Although 
the  authors  in  [6]  incorporate  jump  conditions  in  their  variational  state, 
ment  by  the  addition  of  certain  terms  to  the  classical  functional,  they  do 
not  use  this  modified  form  and  hence  it  is  not  shown  how  effective  it  is. 
Moreover,  since  the  additional  terms  mentioned  above  are  arrived  at  by 
trial  and  error,  their  physical  significance  is  not  immediately  obvious. 

In  this  paper  we  shall  develop  general  variational  statements  for 
harmonic  waves  in  composites,  in  which  the  displacement,  the  stress,  and 
the  strain  in  ono  case,  and  the  displacement  and  the  stress  in  anot  .er,  are 
given  arbitrary  variation.  In  addition  we  shall  develop  all  the  discontinuity 
and  the  quasi- periodicity  conditions  in  a  straightforward  and  systematic 
manner.  Of  course,  we  shall  base  our  developments  on  variational  theorems 
which  have  teen  explored  by  a  number  of  investigators  in  the  mechanics 
literature,  beginning  with  the  work  of  Hellinger  (1914)  [7],  the  unpublished 
thesis  by  Prange  (1916)  [8],  contributions  by  Reissner  (1950-53)  [9,  10], 
works  by  Hu  (1 954)  [  1 1  ],  Washizu  (1  955)  [  1 2] ,  and  a  large  number  of 
other  researchers;  see  the  textbook  by  Washizu  [13]  for  further 
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discussion;  see  also  [14],  In  particular,  regarding  the  jump  conditions 
we  shall  generalize  the  work  by  I-rager  (1967)  [15]  to  include  the  variation 
of  weighted  averaged  tractions  and  displacements  on  a  discontinuity 
boundary.  For  composites  such  a  modification  proves  quite  useful,  since 
the  elastic  constants  of  the  two  constituents  may  differ  substantially  from 
each  other.  Finally,  to  demonstrate  the  effect. vene9s  of  our  results  we 
discuss  in  some  detail,  elastic  waves  propagating  normal  to  the  layers  in 
a  layered  composite.  Our  numerical  results  are  then  compared  with  those 
reported  in  [1,  5,  6]  and  the  exact  solution,  which  comparison  reveals  the 
superiority  of  the  variational  statements  developed  herein;  see  Tables  I 
and  II.  In  fact,  accurate  results  as  those  presented  here,  have  not  been 
obtained  before  by  any  of  the  methods  mentioned  above;  except,  cf  course, 
[4]  where  the  exact  dispersion  relation  is  obtained.  In  addition,  a  second 
example  is  worked  out,  which  illustrates  the  significance  of  the  proposed 
variational  statement  that  includes  a  weighted  averaged  traction  on  interior 
discontinuity  surfaces.  Here,  with  a  suitable  choice  the  weighting 
parameter,  accurate  results  are  obtained;  see  Table  III. 
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2.  ONE- DIMENSIONAL  CASE 


In  order  to  stress  the  essentials  we  begin  by  considering  an 
elastic  medium  whose  properties  vary  periodically  in  the  direction  of  the 
propagation  of  elastic  waves,  i.  e. ,  the  x- direction.  Let  a  be  the 
periodicity- length.  Then  we  have  p  (x  +  a)  =  p  'x),  and  T)  (x  +  a)  =  77  (x), 
where  p  is  the  mass-density,  and  TJ  stands  for  X  +  2p.  when 
dilatation il  waves  are  considered,  and  for  p.  when  shear  waves  are 
considered,  X  and  p  being  the  Lame  coefficients. 

For  harmonic  waves  with  frequency  a),  all  the  field  quantities  are 
proportional  to  e^*^,  where  i  =  \TT7  and  t  measures  time,  and 
therefore  we  have 

H  +  Pll'u  «  «  .  *  =  1*  •  ‘  *  £j  •  (2-D 

where  a  is  the  stress,  u  the  displacement,  and  C  the  strain. 

For  a  periodic  medium  Eq.  (2.  1)  has  periodic  coefficients,  and 
according  to  the  well-known  Floquet  theory,  see  for  example  [16],  it 
admits  a  solution  with  the  property 

g  (x  +  a)  =  g  (x)  eiqa  ,  (2.  2) 

where  g  stands  for  any  of  the  dependent  variables  a,  e  or  u,  and  q 
is  the  overall  wave  number. 
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We  wish  to  establish  general  variational  statements  which  can  be 
employed  together  with  some  appropriate  test  functions  which  may  jr  may 
not  be  continuous  throughout  a  representative  cell,  in  order  to  obtain 
approximate  solutions  to  system  (2.  1).  Since  our  objective  is  to  apply 
these  variational  statements  to  composite  materials  whose  properties  vary 
substantially  from  one  constituent  to  the  next,  we  may  wish  to  use  different 
test  functions  to  represent,  for  example,  the  displacement  in  the  region 
occupied  by  each  constituent  material.  Therefore,  within  a  cell,  the  test 
displacement  function  may  suffer  a  finite  discontinuity  across  the  interface 
between  the  two  materials.  Similar  remarks  apply  to  the  stress  field. 
Hence  a  general  variational  method  must  account  for  such  possible  discon¬ 


tinuities. 

a  a 

Let  us  choose  the  origin  of  the  coordinate  such  that  -  —  s  x  £  — 

defines  a  complete  cell.  In  this  region  let  x  =  x  be  a  point  of  discontin- 

o 

uity  for  a  test  function.  If  the  materi?.!  to  th  *  left  is  completely  disconnected 

from  the  material  to  the  right  of  the  discontinuity  point,  then,  since  we  are 

dealing  with  a  one- dimensional  problem,  we  may  prescribe  arbitrarily  one 

boundary  condition  at  x  =  x  and  one  boundary  condition  at  x  =  x  +  ,  a 

o  o 

total  of  no  more  than  two  conditions;  x  *  =  lim  (x  ±  a),  Ot  >  0.  Hence, 

o  .  o 

a-*  0 

in  general,  one  may  prescribe  arbitrarily  jumps  in  the  value  of  two  field 
quntities  at  x  =  x^  ;  in  the  three-dimensional  case  this  number  increases 
to  six. 


In  addition  to  the  jump  conditions  we  must  'ount  for  the 
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quasi- periodicity  condition  (2.  2)  which  can  be  written  as 


“<!> 


u(--)e-* 


a  (|)  =  (J  (-  |)  e‘qa  .  (2.3) 


As  we  shall  see  below,  these  conditions  can  be  accounted  for  by  either 
the  use  of  a  Lagrangian  multiplier,  or  by  simply  calculating  the  work  of 
the  tractions  at  the  two  boundary  points  of  the  cell. 


General  Variational  Statement  [7-15] 

The  most  general  variational  statement  for  real- valued  field 
quantities  with  discontinuities,  in  which  the  stress,  the  strain,  and  the 
displacement  fields  are  varied  independently,  maybe  stated  as  follows 


I 


a/2 

f 

J 

-?.,c 


pfu  -  a( '  -  f;)]  dx 


or  A  (u  -  U)j- 


a/2 


’-a/ 2 


|  a  <  u) | 

X  =  X 

o 


(2.4) 


where  f  is  the  body  force,  T  is  the  prescribed  value  of  traction  at  the 
end  points,  A  is  the  Lagrangian  multiplier,  and  u  is  the  prescribed 
value  of  the  displacement  at  the  end  points;  *t  each  end-point  we  can 
either  prescribe  the  displacement  or  the  stress,  but  not  both-  Additionally 
in  (2.4)  we  have  used  the  notation 

a  =  ao^+  (1  -  a)  ,  <u)  =  u^-  ,  (2,5) 

where  at  point  xQ,  g^=  g  (x^  ),  and  g^1-  g(xQ),  g  standing  for  either 
<7  or  u;  in  (2,5),  a.  is  an  arbitrary  real  constant,  a  weighting  parameter. 
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whose  significance  will  be  discussed  below.1  Taking  the  first  variation  of 
(2.  4)  we  obtain 

M  =  Ja/2  {[’M-o]  S«  -  [c  -  £]  6a  -  [£  +  pf]6u}& 

I-  1  ia/2 

6  u  -  u  -  U  6 A  j- 

L  J  -a/2 

+  |  <a)  6 u  -  <u)  6a  j-  ,  (2.6) 

X  =  X 

o 

where 

u  =  (1  -  a)  u<3)+  au(1)  •  (2.  7) 

From  (2.  6)  we  observe  that  the  vanishing  of  the  right-hand  side  for 

arbitrary  variation  of  strain,  stress,  and  displacement  in  the  region 
a  a 

-  —  ^  x  ^  —  ,  yields  respectively,  Hooke's  law,  the  definition  of  strain, 

and  the  momentum  equation.  From  the  second  term  on  the  right-hand  side 

of  (2.  6),  moreover,  we  obtain  the  stress  boundary  data  or  the  displacement 

boundary  data,  depending  on  whether  the  traction  or  the  displacement  is 

prescribed  at  a  boundary  point;  in  the  latter  case  the  Lagrangian  multiplier 

A  is  the  traction  corresponding  to  the  prescribed  displacement.  The  last 

term  in  the  right-hand  side  of  (2.  6)  corresponds  to  the  jump  conditions  at 

the  discontinuity  point  x  =  x^  .  At  this  point  one  can  prescribe  arbitrarily 

the  variation  of  no  more  than  two  field  quantities.  To  arrive  at  a  more 

general  case  we  have  implemented  the  weighting  parameter  a  whose  value 

can  be  selected  arbitrarily.  Observe  that,  for  arbitrary  variation 

1  When  two  materials,  1^  ,  0  *  1,2,  meet  at  point  x  ,  the  superscript  is  iden¬ 
tified  with  the  corresponding  material;  i.e.,  a  °is  assigned  to  one  of  the 
constituents. 


of1  OF ,  the  last  term  in  (2.  6)  guarantees  the  continuity  of  the 
displacement  field  across  the  interface  x  -  x  ,  whereas,  for  arbitrary 
variation  of  u  the  stress  is  required  to  be  continuous-  For  a  =  ~  ,  the 
last  term  in  (2.  6)  becomes 

~(6uft)+  6u(l))  (<J*a)-</x))  -  ~i  ( +  5 iP)  (u^.u^)  .  (2.8) 

L*  Cj 

Similarly,  for  Qt  =  0  and  a  --  1  we  obtain,  respectively, 

<c>6uw-  \u>6a'(l),  <a)6uCL>-  <u>6o<3>.  (2.9) 

These  equations  indicate  the  fact  that  at  a  discontinuity  point  different 
quantities  can  be  assigned  arbitrarily.:  Theoretically,  these  equations  are 
equivalent.  In  actual  calculation,  however,  a  proper  choice  for  the  value  of 
g  may  lead  to  more  accurate  results;  we  shall  exemplify  this  later  on. 

Application  to  Composites 

Consider  now  a  composite  consisting  of  layered  elastic  media 

bonded  together,  and  let  us  consider  harmonic  waves  propagating  normal 

to  the  layers.-  For  simplicity  assume  that  a  cell  in  this  composite  consists 

of  two  materials,  ,  P  =  1,2,  where  M1  occupies  the  region 

.  -2  s  x  £  and  £  x  £  ,  and  M3  occupies  the  region 

L  c*  £  +* 

.  —  5  x  £  —  ;  see  Fig.  1,  The  mass  density  and  the  elastic  constantsof  these 
2  2 

1  Note  that  a  is  a  linear  combination  of  and  <j^as  defined  by  i2.  5). 
However,  it.is  a  that  has  independent  variation,  a*:.a  not  C'2 '  and 
o(x )  separately. 
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materials  will  be  denoted  by  p  and  Tf  ,  |3  =  1,  2,  respectively, 

P  P 

where  each  constituent  may  be  inhomogeneous  in  the  direction  of  wave 
propagation. 

To  specialize  the  functional  (2.  4)  to  this  case,  we  observe  that 
this  functional  must  remain  real- valued,  although  the  present  case  the 

field  quantities  are  complex- valued.  Identifying  the  body  force  f  with 

la  .1 

—  (0  u,  're  write 

^  j  ^ 

\  =  I  {  2*»ce*  -  «(«*-  +  P'P'}*1 2 

-a/2 

x=  ±b/2 

(2.  10) 

where  the  superscript  star  denotes  the  complex  conjugate,  and  the  term 
c.  c.  stands  for  the  complex  conjugate  of  the  quantities  which  precede  it. 

The  second  term  in  the  right-hand  side  of  (2.  10)  corresponds  to 
the  constraint  on  the  displacement  at  the  end  points  imposed  by  the  quasi¬ 
periodicity  c  nditior  (2.  3)^  ;  A  is  a  Lagrangian  multiplier.8  Instead  of 
using  the  displacement  constraint,  we  may  replace  this  term  by  a  term 
which  corresponds  to  the  work  of  the  end  tractions.  The  results,  however, 
will  be  the  same,  since  twice  this  work  is  given  by 


1  Here  we  use  the  obvious  fact  that  a  complex- valued  quantity  plus  its1 
complex  conjugate  is  equal  to  twice  of  its  real  part.: 

2  In  the  actual  calculation,  A  must  be  set  equal  to  a  (cj ). 
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°  (!)“*(!) -°(-!)  u*(-!) +  c- 


*  ff(-!)[u*(!)eiqa-  u‘(-l)]  + 


c.  c. 


(2. 11) 


where  (2.  3)  is  used.  The  comparison  of  the  right-hand  side  of  (2.  11) 
with  the  second  term  in  the  right-hand  side  of  (2..  10)  reveals  that  A  =  7 
Taking  the  first  variation  of  Ix  we  now  obtain 
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a/2 


1=  J-a/2 

+  c .  c .  1  dx 

+  {[a(D- A]6u‘(!)  -  [°(-D  -  Ae"‘,a]6u*(-|) 
-  ["*(!)  •u*(-!)eiqi]8A',  +  c  c  } 


+  |  (a)  fiu*  -  <  u)  6a*  +  c.c.:^ 


(2.  12) 


x  =  ±b/Z 


which  for  arbitrariness  of  the  indicated  variations,  gives  all  the  field 
equations  and  the  corresponding  boundary  and  jump  conditions. 

The  functional  1^  may  now  be  specialized  to  yield  other  functionals 
which  may  be  more  suitable  for  calculation  in  a  given  context.;  For 
example,  if  we  assume  that  Eqs.  (2.  1)  are  satisfied  by  the  test 

Cjy  J 

function  for  the  displacement,  1^  reduces  to 
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*.  ■  C  {”£  £  -  ■•*}*  -  *[•*(!)  -  »*(-!Kk,a] 


{("  *)<ui>}  +  cc- 

x  =  ±b/2 


(213) 


whose  first  variation  is 


c 


c. 


} 


dx 


+  ibt  (!)  -  \K(!)  -  |>£  (-1)  -  A.-- ]..*(-!) 

-[•  (!)  - u  (-1) ei,a] 6A* +  '•<•} 

+  {( (”^)>  6“*- <u>  6(^?) +  c  c- }  •  «*•»■ 

x=  ±  b/2 

wV'ere 


<XVa 


duf*> 


(> 


(2.15) 


where  the  superscript  ref<  rs  to  the  corresponding  material,  i.  e.  , 

M1  or  M*  . 

Other  special  cases  are  obvious.  For  example,  if  the  test  function 
is  chosen  such  that  the  quasi- periodicity  condition  is  automatically 
satisfied,  the  second  term  in  the  right-hand  side  of  (2.  13)  drops  out.  If. 
moreover,  the  test  function  is  continuous  throughout  the  cell,  the  last 
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term  in  the  right-hand  side  of  (2. 13)  vanishes.  In  this  very  special  case 
we  then  arrive  at  the  usual  Rayleigh  quotient  which  defines  the  frequency. 
This  is  the  case  which  was  exemplified  in  detail  in  [5,  6].  Observe  that 
the  corresponding  approximate  expression  for  the  stress  field  obtained  in 
this  manner  would  suffer  a  discontinuity  at  the  interface  of  the  two  constitu¬ 
ent  materials.  The  results  are,  therefore,  very  poor;  see  Table  II.  If, 
on  the  other  hand,  different  test  functions  are  used  to  represent  the  displace¬ 
ment  field  within  the  region  occupied  by  each  constituent  material,  the  last 
term  in  (2.  13)  must  be  retained.  This  may  then  result  in  a  more  reasonable 
approximate  expression  for  both  the  displacement  and  the  stress  fields. 

We  therefore  wish  to  stress  here  that,  if  one  wishes  to  vary  only  the  dis¬ 
placement  field,  as  is  done  in  [5,  6],  one  then  must  use  this  latter  approach, 
rather  than  the  usual  Rayleigh  quotient. 

In  most  calculations  it  is  much  easier  to  use  the  same  set  of  test 
functions  throughout  the  entire  cell.  To  arrive  at  sufficiently  accurate 
results,  another  variational  statement  must  be  employed,  in  which  both 
the  displacement  and  the  stress  fields  can  be  directly  and  independently 
approximated.  In  the  following  we  shall  discuss  this  case. 

Another  General  Variational  Statement  [7-151 

Let  us  denote  by  D  the  inverse  of  the  elasticity  coefficient  t) ,  i.  e.  , 

D  -  ~  •  (2.16) 

The  most  general  variational  statement  for  real-valued  field  quantities 
with  discontinuities,  in  which  the  displacement  and  the  stress  fields 
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(but  not  the  strain  field)  are  varied  independently,  may  be  stated  as 
follows;  compare  with  Eq.  (2.4): 


ra/2  r  i 

J  =  —  Da3  -  p  f  u  - 

-a/2  LZ 

du  "1 

dx  +  | 

Tu, 

,  n102 
or  A(u-O)} 

-  a/2 

+  1  a  <  u )  \ 

(2.  17) 

v  / 

x=x 

o 

where  all  the  quantities  are 

defined  as  before. 

The  first  variation  of 

(2.  17)  is 

J  6a  + 

[£♦ 

pto2 

u  J  6  u  ^  dx 

-  {[« -  t]6°  ■ 

or  £a  -  A 

|  6u 

_  ^  a/2 

r  u  -  uIsa} 

L  J  -  a/2 

-  |  (a)  fiu  -  (u) 

6  a  | 

(2.  18) 

x-x 

o 


which  yields,  for  arbitrary  variation  of  the  indicated  quantities,  all  the 
field  equations,  and  the  boundary  and  jump  conditions.: 

Observe  that  if  continuous  displacement  and  stress  test  functions 
are  used  in  (2.  17),  the  last  term  drops  out-.  Nevertheless,  the  corres¬ 
ponding  variational  statement  will  yield  reasonable  approximate  expressions 
for  the  displacement  and  the  stress  fields. 

To  apply  (2.  17)  to  elastic  waves  in  composites,  we  proceed  as 
before,  and  consider  the  functional 
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J  =  i  —  D OO  +  —  pa)  uu  -  ff  — —  +  c.  c.  )  dx 

l  J  L  2  2  dx 

-  a/2 


ra/2  f  j 


.2 


du^ 


■} 


{A[u*(l)  -  u*(-l)e'”,i]} +  ia  <u”'}  *cc- 

x=±b/2 


(2-19) 


whose  first  variation  is 


C  {[»-£]•••*  [2 
-  {[*(!)-  A]6u’(!)- [•’(-!)- 

-  [u  (1)  - u  (-!)e",a]6A* +  c  c- 

-  •f<o)6u*-  (a)  6a*  +  c.c.  1 

^  x=±  b/2 


* 


+ 


c.  c. 


} 


dx 


} 


(2.  20) 


It  is  clear  that  the  vanishing  of  6^  for  arbitrary  variation  of  the  indicated 
quantities,  guarantees  the  satisfaction  of  the  field  equations,  the  quasi¬ 
periodicity  conditions,  and  the  continuity  of  the  displacement  and  the  stress 
across  the  interface  of  two  materials  within  the  cell.. 

Let  us  now  consider  numerical  examples  before  discussing  the 
general  three-dimensional  case-. 

Illustrative  Examples 

I.  As  our  first  example  let  us  illustrate  the  effectiveness  of  our 
variational  theorem  corresponding  to  functional  (2..  19).  We  assume  that 
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each  constituent  in  the  composite  is  homogeneous  and  isotropic. 


As  our  test  functions  let  us  consider 


u  = 


±N 

£ 

ft=0 


U  e 
n 


i(Q+  2irn)£ 


±N 

E  S  e 
n=o  n 


i(Q+  2irn)£ 


(2.  21) 


where  Q  =  qa  ,  £  =  —  ,  and  U 

a  n 

which  are  to  be  calculated;  (2.  21)^ 


and  S  are  the  Fourier  coefficients 
n 

is  the  test  function  used  in  [5,  6]. 


Since  the  above  test  functions  are  continuous  throughout  the  cell, 
and  moreover  satisfy  the  quasi- periodicity  conditions,  only  the  integral 
in  (2.  19)  survives.  Substituting  from  (2.  21)  into  this  integral  and  carrying 
out  the  indicated  integrations,  we  arrive  at 


J  (N) 

l 


±N 

E  Apoj3  U  U 
„  “f  n  L  n  i 

n,  m=U 
n^m 


m 


ADS  S 
n  m 


] 


sin  it  (n 


m)  a 


ir  (n  -  m) 


±N  _ 

+  E  |PW2u  U*  +  DS  S*  +  -(Q+  2vn)(S  U*  -  S*  U  )]  ,  (2.22) 

L  nn  nna  nnn  n  J 

n=U 


where  the  following  notation  is  employed: 


Ap  =  p  -  p  ,  AD  =  D  -  D 

3  1  3  1 


p  =  n  p  +  n  p  , 

1  1  3  3 


D  =  "l  D1  +  n3  Da  * 


(2.23) 


a  -  b 

n  =  - 

i  a 


b 

n  =  — 
s  a 


To  begin  with  we  immediately  observe  that,  for  N  =  0,  (2.  22) 


yields 
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(d  =  Q/fpD)^3 


(2.  24) 


which  corresponds  to  the  non- dispersive  results  that  can  be  obtained  by 
the  usual  method  of  calculating  the  effective  mass-density  and  elastic 
modulus. 

For  N  a  1,  we  introduce  the  notation 


U 

S 


— 

1  u  „  , 

U  M  ,  , 

.  . .  ,  U  ,  . .  . , 

l  -N  ’ 

-N+l 

o 

f  s 

s  ,  , 

c 

l  -N  ’ 

-N+l  ' 

o’  ■■■’ 

(2.  25) 


where  superposed  T  denotes  the  transpose,  and  write  (2.  22)  as 


where 


H  =  diag  ~  |  Q  -  2irN  ,  Q  +  2tt  <n  +  1)  ,  ,  .  .  ,  Q,  .  ,  .  ,  Q  +  2irN  j 


(2.  27) 
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and  whe re  the  matrix  ♦  is  obtained  by  replacing  aJ2  p  in  the  matrix 
Q  by  D,  and  a.’3  Ap  by  AD,  respectively.  From  che  stationary  con¬ 
dition  imposed  on  (2.  26)  we  now  immediate  arrive  at 

fl  U  +  H  S  =  0  ,  H*U  +  *S:--0  , 

JV*  JNJ  ^  <“N«*  j  J«S/ 

which  yield 

S  =  ,  r  H*  -  H  1  0 1  U  =  0  (2.  28) 

AJ  aj  <V  I  AJ  AJ  AJ  JV  I  AJ  <V 

The  second  equation  in  (2.  28)  provides  the  eigenfrequencies  and  the 
corresponding  eigenvectors,  and  the  first  equation  then  gives  the  Fourier 
coefficients  for  the  stress  field.  The  characteristic  equation  therefore  is 

det  f  H*  -  *  H"1  ol  =  0  (2.29) 

L  ~  ~  ~  ~  J 

For  a  given  value  of  the  wave  number  Q,  the  roots  of  this  equation  give 
the  corresponding  frequencies. 

Table  I  gives  typical  results  with  the  corresponding  exact  values  for 

1  Pa  \ 

n  =  n  =  —  ,  —  =  3,  and  the  indicated  values  of  —  .  In  this  table  the 

i  a  2  p  T) 

l  'i 

values  of  the  dimensionless  frequency, 

V  =  aw(p/T))^a,  T)  =  nx\  +  n3  >  (2.30) 

are  listed  for  the  indicated  values  of  Q  and  N,  and  for  various  eigen- 
modes.  As  is  seen  from  this  table,  our  crudest  approximation  which 
still  manifests  dispersion  and  which  corresponds  to  N  =  1,  (i.  e.  ,  3  plane 
waves),  gives  for  the  first  mode,  results  which  are  extremely  accurate. 
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For  higher  modes  N  must  be  taken  greater  than  1,  although  even  for 
N  =  1,  the  secoixi  and  third  modes  obtained  from  (2.  29)  are  still  quite 
accurate;  we  shall  discuss  this  below.  As  is  seen  from  Table  I,  for 
N  2  2,  again  extremely  accurate  results  are  obtained.  Indeed,  to  our 
knowledge,  no  other  approximate  method  has  as  yet  produced  such  superior 
results  with  so  little  computational  effort. 

To  demonstrate  this  last  assertion,  we  have  compared  some  of 
our  results  with  those  corresponding  to  the  methods  reported  in  [l]  and 
in  [5,  6],  in  Table  Ila.  Observe  that,  even  for  N  =  5,  (i.  e.  ,  11  plane 
waves),  the  method  of  Kohn,  et  al.  [6]  gives  poorer  results  than  those 
deduced  from  Eq.  (2.  29)  for  N  =  1,  (i.  e.  ,  only  3  plane  waves/,  while  for 
a  fixed  N,  Eq.  (2.  29)  requires  less  computational  effort  than  the  corres¬ 
ponding  equation  in  [6];  this  is  b^ause,  in  Eq.  (2.  29)  matrix  H  is 
diagonal,  V/hereas  in  [6]  the  characteristic  equation  has  the  form 
det  |  A  -  t/a  B  |  =0,  in  which  neither  A  nor  B  is  diagonal. 

In  Table  lib  the  values  of  the  frequency  parameter  v  for  higher 
modes,  as  given  by  Eq.  (2.  29),  are  compared  with  those  obtained  in  [6]. 

We  note  that  for  N  =  2,  for  example,  the  characteristic  Eq.  (2.  29) 
is  a  5x5  determinant  which  can  yield  only  the  first  five  eigen- 
frequencies.  As  a  rule,  only  the  first  eigenfrequency  is  reasonably 
accurate.  A  very  remarkable  feature  of  the  present  method,  howover,  is 
that,  even  the  higher  frequencies  obtained  in  this  manner,  are  quite 
accurate,  whereas  this  is  not  the  case  .or  the  method  proposed  by  Kohn, 
et  al.  [6].  Table  lie  illustrates  this  fact.  In  this  table  all  eigenfrequencies 
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are  reported  for  N  =  2. 


EL.  For  our  second  example  we  shall  illustrate  the  variational  state¬ 
ment  (2.  13)  and  especially  the  significance  of  the  weighting  parameter  a. 
Again  we  consider  waves  propagating  normal  to  the  layers  in  a  layered 
medium.  Here  we  shall  use  a  different  displacement  field  in  each  con¬ 
stituent,  and  set 

(i)  „  -  '  a  a  ^  _  b 


±  N  i  (Q  +  2im)— 

u-  •  =  L  U  e  a 

n=0  n 


-  2  *  X  S 


u 


(3) 


±N 


m=0 


V  e 
m 


i  2  it  m  ~ 
b 


j  b  ,  a 

and  —  s  x  ^  — 

2  2 


b  <  <  b 

..  —  s  X  S  — 

2  2 


(2.31) 


Whil^  u^  ^  satisfies  the  quasi- periodicity  condition,  u^  does  not,  and 
therefore  (2.31)  represents  a  crude  approximation.  Nevertheless,  if  the 
weighting  parameter  a  is  chosen  appropriately,  reasonable  results  can 
be  obtained.  We  shall  not  report  the  detailed  algebra  here,  and  only  give 


the  final  equation: 
±N 


*  n  r  . 

~  y  .S 
1  "n  y 

n.m=0  [_ 


n,m= 

nj^m 


(Q+  2irm)  (Q+  2Trn)  +  +  -==-  ira(m  -  n) 

0  y 


3  A  sin  it  (m  -  n)— 

V  .  4  a ,  .3  a  TT  * 


3 


/  v  u 

ir  (m  -  n)  m  n 


±N  (f  i  V*nT\  r.,  .3  |/an  01  ) 

£  ]  (Q+  2irn)an  .  —2  U*U  +  -  — a-  V*V  } 

n=o!Ly  1  9  J  n  n  Lnay  0  J  n  n) 


±N  ,r, 

+  2  £  Jr  & 

n1m=0  y  L  " 


2-rrnqy 


+  (1  -  a)  (Q+  2'.rm)J  sin|  ir(n-  m-^)  -  Q 


l  IT 


_b 

2a 


I- 


+  U  V 
m 


(2.  32) 
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« IJtKIMWUWJ 


Off  -W 


■rj^v  -wn'vx*  “*»T' 


where  y  =  fj/tj,0  =  p/p,y  =  n+  yi  ,  and  0  =  n  +  0n  . 

3  1  3  1  1  12 

Table  Ilia  gives  the  frequency  parameter  V  for  0  =  3,  nx  =  n^  =  ■—  , 

N  =  5,  and  indicated  values  of  Q  and  a.  As  is  seen,  a  can  have  a 

significant  effect  on  the  accuracy  of  the  result.  It  should  be  remarked 

that,  in  view  of  the  discontinuity  conditions  which  are  involved  now,  the 

algebra  is  more  cumbersome  than  in  the  previous  example.  Another  point 

that  should  be  stressed  is  that  a  =  —  is  not  necessarily  a  good  choice  for 

this  parameter,  as  is  seen  from  Table  III;  for  a  =  ~  our  general 

c* 

variational  statement  (2.  13)  reduces  to  the  case  considered  in  [6],  The 
optimum  value  of  a  occurs  at  a  =  2,  and  it  appears  that  this  optimum 
a  is  independent  of  the  value  of  y  =  t)  /f ^  ,  for  fixed  values  of  nx 

and  0,  and  for  N  ■  5.  Further  calculation  revealed  that  for  sufficiently 
lar*je  N,  say,  N  *  3,  the  optimum  value  of  or  is  always  2,  independently  of 

the  values  of  the  other  parameters.  For  small  N,  say,  N  <  2,  on  the  other 

hand,  optimum  a  changes  with  the  other  parameters.  This  is  exemplified  in 
Fig.  2,  where  tl>*  variation  of  v  is  plotted  against  a  for  indicated  values 
of  the  other  parameters.  Note  that  for  N  ■  0,  the  optimum  value  of  a  is  1.0, 
whereas  for  N  ■  2,  the  optirosn  value  Is  2.1.  It  is  remarkable  that  even  for 
N  •  0  (i.e.,  u(l>  *  UQe^  and  u^  ■  , ,  a  m  1.0  gives  v  *  0.280  as  com¬ 
pared  with  the  exact  value  v  ■  0.271.  For  N  ■  2  and  a  *  2.1,  on  the  other 

hand,  v  •  0.274. 
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3.  THREE-DIMENSIONAL  CASE 


We  shall  first  develop  the  relevant  general  variational  statements 
and  then  specialize  the  results  for  application  to  harmonic  waves  in 
composites  with  periodic  structure.  As  before,  we  shall  be  exclusively- 
concerned  with  the  linear  elasticity  theory. 

Using  a  rectangular  Cartesian  coordinate  system  we  denote  the  dis¬ 
placement,  the  strain,  and  the  stress  components  by  u.  ,  €  ,  and  a  , 

J  J  J 

respectively,  and  assume  that  the  body  forces  f.  measured  per  unit  mass 

are  given  throughout  the  volume  V  occupied  by  an  elastic  body  with  mass- 

density  p  and  elasticity  coefficients  C„  .  (x),  where  x  is  the  position 

jk-tm  ~  ~ 

vector  of  a  typical  point  in  V;  j,  k,  t,  m=l,2,3.  In  the  general  anisotropic 

case  C  is  symmetric  with  respect  to  the  exchange  of  j  am1  k, 

jk-om 

l  and  m,  and  jk  and  -tm.  For  the  isotropic  case,  moreover,  we  have 

C„  ,  =  X6..  6,  +  2u  6..  6,  ,  where  X  and  n  are  the  Lame  coeffic- 

jk-tm  jk  Im  r  jl  km 

ients  (they  may  depend  on  x),  and  6  is  the  Kronecker  delta;  we  need 
not  confine  our  discussion  to  the  special  isotropic  case. 

The  field  equations  are: 


<*...+  Pf.  =  0 

jk,  k  j 


fjk  "  U(j,  k) 


CTjk  ^jktm 


U(j,  k) 


-r  (u  +  u.  .) 
i.  i,  k  lc,j' 


(3.  1) 


where  a  comma  followed  by  an  index  letter  indicates  differentiation  with 
respect  to  the  corresponding  coordinate  variable,  and  repeated  subscripts 


2. 


are  to  be  summed  from  1  to  3. 


To  characterize  the  boundary  'ta  wc  denote  by  n.  the  exterior 

unit  normal  on  the  regular  surface  o  which  bounds  the  body,  and  assume 

that  at  a  typical  point  on  S  certain  components  of  the  displacement  vector, 

together  with  the  complementary  components  of  the  tr  actio/*  vector,  are 

prescribed.  We  denote  by  T.  the  prescribed  traction  components, 

1 

and  by  the  prescribed  displacement  components; 

* 

henceforth,  singly  and  doubly  underlined  subsc-^pt  letters  refer,  respec¬ 
tively  to  the  prescribed  tractir  and  displacement  components  at  a  point 
on  S.  We  thus  have 

O’.,  n.  =  T,  ,  u.  =  U.  ,  on  S  (3.2) 

ik  K  1  11 

In  our  general  variationa1  statements  we  shall  admit  test  functions 
which  may  have  finite  discontinuities  across  a  finite  number  of  isolated 
surfaces  within  the  volume  V.  We  shall  denote  the  collection  of  these  dis¬ 
continuity  surfaces  by  S,  and  obterve  that,  at  a  point  on  X,  three  boundary 
data  can  be  prescribed  for  the  material  on  one  side  of  the  surface,  and 
three  boundary  data  for  the  material  on  the  "ther  side,  a  total  of  no  more 
than  six  conditions.  Hence,  a  jump  in  the  value  of  at  most  six  quantities 
may  be  prescribed  arbitrarily  at  a  point  on  S- 

With  the  above  preliminaries  out  of  the  way  we  now  consider  a 
variational  statement  in  which  the  displacement,  the  strain,  and  the  stress 
can  be  given  arbitrary  variation,  as  follows..  Consider  the  functional 
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K 


^  (  2  Cjktm  €jk  €lm  "  P  f j  U j  "  ffjk  L  *jk  '  “(j,  k)  J}  dV 


I 


T.  u.  dS 

S  1  2 


f  T.  f  u.  -  U.  IdS  -  f  T.  <u.)  dE  , 


(3.3) 


where 


T.  =  aTfe)-  (1  -  a)  t[x)  ,  <u.)  =  ufe)-  u.(l?  on  E  .  (3.4) 

J  J  J  J  J  J 


Here  the  superscripts  are  defined  in  the  following  manner:  The  discontin¬ 
uity  surface  E  divides  the  body  into  subregions  with  at  least  a  part  of 
E  as  a  common  surface.  Consider  two  such  adjacent  subregions,  say, 
regions  1  and  2,  and  the  field  quantity  g<  At  a  typical  point  on  E,  the 
quantity  g  ,  p  1,2,  then  denotes  the  limiting  value  of  g  as  the  con¬ 
sidered  point  is  approached  from  the  interior  of  subregion*  p. 

In  Eq.  (3 -  4; ^ ,  a  =  C*(x)  is  a  real- valued  function  defined  on  E;  this 
function  can  be  chosen  arbitrarily  to  expedite  the  numerical  calculation  in 

the  given  context.  [For  composites,  a  is  assigned  to  one  of  the  constituents. 3 
If  we  now  observe  that  the  exterior  unit  normal  to,  say,  subregion  1, 
is  the  negative  of  the  exterior  unit  normal  to  subregion  2  at  a  typical  point 

on  E,  the  first  variation  of  (3.3)  becomes 

5K  ivlLCjktm 'lm  ‘  CTjk]  fiCjk  *  L'jk  "  “(j,  k)J8°jk 

’  [°jk.k  +  Pfj]6“j}dV 

+  Js["iknk-  TJvs-  Uv  uJ6ys 

*  For  composites,  8  refers  to  the  corresponding  material,  iP . 
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(3.5) 


where 


u. 

J 


(1  -  a)  u<a)+  Ofufl)  on  S 
J  J 


(3.6) 


We  observe  that  for  arbitrary  variation  of  strain,  stress,  and  displace¬ 
ment  in  V,  we  obtain  Hooke's  law,  the  definition  of  strain,  and  the 
momentum  equations,  respectively.  On  the  boundary  S,  moreover, 
arbitrary  variation  of  u.  and  T.  gives,  respectively,  the  prescribed 

2m  *1 

traction  and  displacement  boundary  data.  The  last  term  in  (3.  5)  corres¬ 
ponds  to  *he  jump  conditions  on  the  discontinuity  surface  £•  It  indicates 
that,  if  weighted  averaged  displacement  and  tractions  defined  by  Eqs.  (3.  6) 
and(3.4)j  ,  respectively,  are  given  arbitrary  variation,  then  the  continuity 
of  the  traccion  and  displacement  components  is  guaranteed.  Note  that, 
similarly  to  the  one- dimensional  case,  special  forms  of  these  weighted 
averaged  quantities  can  be  obtained  by  the  suitable  choice  for  the  arbitrary 
function  ol. 

Before  proceeding  to  apply  the  variational  statement  (3,  3)  to 
harmonic  waves  in  composites,  let  us  state  another  general  variational 
theorem  in  which  only  the  stress  and  displacement  are  given  arbitrary 
variation.  To  this  end  consider  the  functional 


Hi 


where  B.  C. 


D.,  ,  a.,  a,  + 

jk-om  jk  Im 


pf.  u.  -  a.,  u...  ,  ,  \  dV  +  B.  C. 
J  J  Jk  (j,k)J 


+  J.C. 


(3.7) 


and  ,.T .  C.  denote,  respectively,  the  boundary  and  the  jump 
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conditions  as  given  by  the  last  two  terms  (with  minus  sign  changed  to  plus) 


in  the  right-hand  side  of  (3.3),  and  D  is  the  elastic  compliance 

matrix  obtained  by  inverting  the  matrix  C  .  The  first  variation  of 

jk-tm 

L  is 


6L  =  Jvt[Di 


jktm  f-m 


"(j,  k)]  6°jk  +  [°jk.  k  *  p  Jj]  6uj  }  dV 


+  B.C.  +  J.C.  ,  (3.8) 

where,  again,  B.  C.  and  J.  C.  correspond  to  the  negative  of  the  last 
three  terms  in  (3.  5). 


Application  to  Composites 

For  the  sake  of  simplicity  in  presentation,  let  us  consider  a 
composite  consisting  of  a  densely  packed  and  completely  bounded  collection 
of  unit  cells  in  the  form  of  parallelpipeds;  our  results,  however,  apply  to 
composites  with  different  periodic  structures,  without  any  additional 
difficulties.  Consider  a  representative  cell,  and  let  its  three  edges  be 
defined  by  three  vectors  $ ,  p  =  1,  2,  3-  We  shall  assume  that  the  cell 
with  volume  V  and  surface  S  consists  of  different  material  constituents 
which  are  separated  from  each  other  by  the  interior  surfaces  £ «  We 
observe  that,  because  of  the  periodic  structure,  we  have  g  +  tP'j  =  g(x), 
where  g  stands  for  any  of  the  material  properties. 

We  shall  consider  harmonic  waves  of  fre  juency  oo  and  vave 
vector  propagating  in  this  composite..  Because  of  the  periodic  struc¬ 

ture  of  the  medium,  the  field  equations  have  periodic  coefficients,  and 


therefore,  according  to  the  Floquet  theory,  they  admit  solutions  of  the 
form 

f  (*  +  exp{iq.t?}  ,  (3.9) 

where  f  stands  for  any  of  the  field  quantities-.  In  particular,  for  x  on 
the  surface  S,  Eq.  (3.9)  defines  the  quasi- periodicity  conditions  for  the 
displacement  or  tractions,  depending  on  which  quantity  is  identified 
with  f.  For  the  sake  of  simplicity  in  presentation  let  us  denote  three 
faces  of  the  parallelpiped  which  intersect  at  a  common  point,  by  S'  ,  and 
designate  a  typical  point  on  S'  by  £  .  Then  it  is  clear  that,  to  each 
point  £  on  S'  there  corresponds,  for  a  suitable  tP  ,  a  point  £  +  lP 
on  S  -  S' .  The  quasi-periodicity  condition  then  becomes 

“jU  +  ip)  *  «j<£>  • 

(3.  10) 

TjU  +  ip)  =  -Tj<£>  “p^Vk}  •  £  on  s'  ■ 

where  the  minus  sign  in  the  right-hand  side  of  the  last  equation  occurs, 
since  the  exterior  unit  normals  on  S  at  £  and  £  +  tP  ,  are  oppositely 
oriented 

As  in  the  one-dimensional  case,  we  shall  employ  the  first  equation 
in  (3.  10)  as  constraints  on  the  displacement  field. 

To  arrive  at  a  variational  statement  in  which  the  displacement, 
the  strain,  and  the  stress  are  varied  independently,  we  now  identify  the 
body  forces  by  —  go  8  u-  ,  and  from  (3.3)  obtain 

u  J 
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K  =  f  Uc..f  C..  €*  -  ~P0)au.  uf  -  a.,  f  t*  -  u*  ..1+  c.c.\  dV 

x  L  2  jktm  jk  Im  2  j  j  jk  L  jk  (j,  k)  J  J 


-  J  '  {a^  [  «*  (£  +  -  «*(£)  exP  (  •  i(*k  ^k  *]  +  c-  c>  }  dS 

s 

-  J  ^  T  j  (  u  *  )  +  c .  c .  j-  d  S 


(3.11) 


where  all  the  terms  are  as  defined  before.  Now,  taking  the  first 
variation  of  Kx  we  obtain 


8KrJ 

y  iLCjkl.m  ajkJ  ^  *jk  “ 

Lv 

-  u(j.k)]6V 

-  [  a..  .  +  pto3  u.  |  6  u  *  + 

C.  C-: 

\  dV 

L  jk,  k  r  j  J  j 

J 

+ 

ft 

fs,{[Tj(i  +  ip)-  A j ] 6 “ * 

U* 

+  [  Tj(^)  +  Aj  exP  f 

6u* <£) 

-  [uj(i  +  -  V£>  exP  1 1\ 

0  6Aj*  +  C  C‘  }  dS’ 

+ 

ft 

[  {  <T.>  6  u  *  -  <u.  >  6T*  + 
k  J  J  J  3 

c.  c.  j- 

dE  (3,12) 

As  is  seen,  the  vanishing  of  the  first  integral  in  the  right-hand  side  of 
(3.  12)  for  arbitrary  variation  of  the  indicated  field  quantities,  yields  all 
the  field  equations.  The  vanishing  of  the  second  integral  shows  that  the 
Lagrangian  multiplie r  A.  is  given  by 

Aj  =  V£  +  £P>  =  -Tj<l>  e*p  £  tqk  4.P}  , 


?7 


Thus  the  second  term  in  the  right-hand  side  of  (3.  11)  is  twice  the  potential 
of  the  surface  tractions  which  are  constrained  to  satisfy  the  quasi- period¬ 
icity  condition  (3.10)  .  Finally,  the  vanishing  of  the  last  integral  in  the 

right-hand  side  of  (3.  12),  for  arbitrary  variation  of  the  weighted  averaged 
displacement  and  tractions,  as  defined  on  £  by  Eqs.  (3.4)^  and  (3.  6), 
guarantees  the  continuity  of  the  tractions  and  the  displacements  on  the 
discontinuity  surface. 

Functional  (3.  11)  can  now  be  specialized  in  the  same  manner  as 
discussed  in  connection  with  the  one- dimensional  case.  For  instance,  if 
we  wish  to  vary  only  the  displacement  field,  we  assume  that  the  last  four 
expressions  in  (3.  1)  are  satisfied,  and  modify  accordingly  the  integrand 
of  the  first  integral  in  the  right-hand  side  of  (3.  11). 

In  a  similar  manner  we  can  immediately  write  down  the  functional 
from  the  functional  L,  which  yields,  for  harmonic  waves  in  com¬ 
posites  with  periodic  structure,  the  appropriate  variational  theorem  in 
which  the  displacement  and  the  stress  fields  are  given  independent 
variation.  In  this  manner  we  obtain 


U-K 


>jt  ,  1  9  #  5j( 

»  0.,  o ,  +  -r  pur  u.  u.  -  ,o  u.  , .  + 

jktm  ]k  to  2  ]  J  jk  (j,  k) 


:.c.} 


Q.  C. 


(3. 13) 


where  Q.  C.  and  J.C.  stand  for  the  quasi- periodic  and  the  jump  con¬ 
ditions,  respectively,  which  are  the  negatives  of  those  occurring  in  the 
expression  for  .  As  can  be  immediately  verified,  the  vanishing  of  the 


first  variation  of  the  L  ,  for  arbitrary  variation  of  the  appropriate 
quantities,  gives  all  the  field  equations,  the  quasi- periodicity  conditions, 
and  the  jump  conditions- 

Application  to  Fiber- Reinforced  Composites 

As  our  final  illustrative  example,  we  consider  harmonic  waves 
propagating  normal  to  the  fibers  in  a  fiber- reinforced  elastic  composite 
whose  fibers  are  approximately  rectangular  in  cross  section,  and  are 
placed  parallel  to  each  other  in  a  periodic  manner.  The  cross  section  of 
a  representative  cell  is  shown  in  Fig.  3.  We  consider  waves  propagating 
in  the  x1 -direction. 

The  general  variational  theorems,  Eqs.  (3.4)  and  (3.  7),  maybe 
applied  directly  here.  As  an  illustration  we  shall  apply  (3.  7)  which 
permits  variations  in  the  stress  and  displacement  fields.  The  other  case 
can  be  handled  in  a  similar  manner. 

With  the  geometry  defined  in  Fig.  3  we  therefore  write 


a  /2 

r 1 

a  n 

r 3 
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a 
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-a  /2 
a 


_  _  _  iqa 

+  ^(t'  x.)[’,»*('T-  \)‘  u2*(-f  ’  O*  l]  +  c  c  }' 


+  I.  ;,Ha(v  *f)  <ua*(v  ±f)>  +  °„(v  *‘#)<u>*(v  *T> 

-  b  /  2 

+  c .  c .  ^  dx^ 


b  n 

r  3  f 

-  /  bi 

f  { 

a  l±  —  »  x 

-b  /2  ^ 

ii  \  2 

2 

+  c.  c.  }•  dx 

J  2 

P.  r.  6 

,  T)  =  1,  2 

(3.  14) 


where  is  the  elastic  compliance  matrix  for  the  plane  strain  case, 

which  for  isotropic  materials  becomes  D_  ,  =  -^-^6  6  -  — :  6.  6.  1, 

y  Py6rj  2pL  p6  y?)  2(X  +  n)  py  fitjJ’ 

where  6-  is  the  two-dimensional  Kronecker  delta. 

0y 

The  second  integral  in  functional  (3.  14)  corresponds  to  the 
periodicity  conditions 


“aU'  f )  *  us(v  -f)  1  p  =  *•  2  • 


(3.15) 


whereas  the  third  integral  pertains  to  the  quasi- peridocity  conditions 

ua(T'  %  )  =  ua(-  ~t  ’  *»)  e‘qai  :  <3’16) 
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in  these  integrals  the  Lagrangian  multipliers  are  identified  with  the 

P 

corresponding  '  rations.  The  last  two  integrals,  moreover,  correspond  to 
jump  conditions,  if  different  test  functions  are  in  each  region  for  the 

displacement  and  stress  fields.  As  we  have  pointed  out  before,  the 
inclusion  of  these  jumps,  although  it  may  improve  on  the  detailed  local 
variation  of  the  field  quantities,  results  in  a  more  complic'ted  numerical 
procedure.:  For  this  reason  it  may  be  desirable  to  use  continuous  test 
functions  for  the  displacement  and  the  stress  components,  in  which  case  the 
last  two  integrals  drop  out.  Simple  test- functions  of  this  kind  are 

up(xi  ’  xa)  =  £  Up^m)  exP  {*  [<Q  +  2lrn)£x  +  2irm£3]}  « 


n,  m=0 


(3.  17) 


V(Xi’Xa)=  ^  S^)  exp  {  i  [  (Q  +  2tm)  ^  +  2nrn4a  ]}  , 

n,  m=0 


where 


Q  =  i3,  •  «i  =  T  •  { 


3  a 


Observe  that  test- functions  (3.  17)  satisfy  the  continuity  conditions,  the 
periodicity  and  the  quasi- periodicity  requirements  defined,  respectively, 

by  (3.  15)  and  (3.  16).  These  functions,  however,  do  not  permit  for  the  dis- 

b 

1 

continuity  of  the  stress-components  O’  on  the  interfaces  at  x  =  ±  —  ,  and 

b  32  l  2 

a  on  the  interfaces  at  x  =  ±  —  .  Since  no  exact  solution  exists  for  this 
11  a  2 

problem,  it  is  not  possible  to  easily  assess  the  accuracy  of  the  results.  For 
this  reason  we  postpone  further  discussions  until  an  assessment  of  their 


accuracy  is  made  by  means  of  a  mure  refined  procedure. 
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TABLE  I 

Frequency  Parameter  v  =  au)  ( p/  TJ  )^a  for  Indicated  Values 

%  b  1  P 

of  Q  =  qa,  N,  and  —  ,  and  for  —  =  —  ,  and  —  =  3. 


FIRST  MODE 


N  =  1 


0.3987 
0.  7893 
1.  4925 
1.  8744 


SECOND  MODE 


N  =  2 


5.  786 
5.  194 
4.  506 
4.  122 


EXACT 


5.  795 
5.  205 
4.  516 
4.  125 


T)  /r)  = 
'a  'i 

50 

N  =  2 

EXACT 

1.954 

1.958 

1.902 

1.907 

1.791 

1.797 

1.722 

1.725 

FIFTH  MODE 


N  =  5 


12.  51 

12.  98 

13.  66 

14.  05 


EXACT 


12.  53 

13.  00 

13.  68 

14.  10 


50 

N  =  5 

EXACT 

6.  45 

6.  49 

6.  54 

6.  58 

6.  76 

6.  78 

7.  00 

7.  03 

TABLE  Qa 

Frequency  Parameter  V  -  &U>  (p/ 1})^2  as  given  by  Eq.  (2.29),  by  Kohn, 

Krumhansl,  and  Lee  [6]  ,  and  by  Sun,  Achenbach,  and  Herrmann  [l]  , 

b  1  Pa 

for  Indicated  Values  of  Q  =  qa  and  N,  and  for  —  =  —  ,  —  =  3,  and 


FIRST  MODE 

Q 

EXACT 

Present  Method 
Eq.  (2.29) 

Kohn,  Krumhansl, 
and  Lee,  Ref.  [6] 

Sun,  Achenbach 
and  Herrmann 
Ref.  [1] 

N  =  1 

N  =  1 

N  =  5 

0.5 

0.  098 

0.  098 

0.  23  7 

0.  107 

0.  099 

1.  0 

0.  193 

0.  193 

0.  460 

0.  210 

0.  196 

2.0 

0.354 

0.354 

0.  854 

0.3  83 

0.381 

3.0 

0.  434 

0.430 

1.  083 

0.467 

0.  546 

4.  0 

0.388 

0.3  80 

1.  034 

0.  421 

0.687 

5.  0 

■ 

0.  244  | 

0.  238 

0.  694 

0.  267 

0.  805 

6.  0 

0.  056 

0.  055 

0.  168 

0.  062 

0.  901 

TABLE  nb 


Frequency  Parameter  V  =  a<o(p/l))^3  for  Higher  Modes,  as 
given  by  Eq.  (2.  29),  and  by  Kohn,  Krumhansl,  and  Lee  [6]  , 
for  Indicated  Values  of  Q  =  qa  and  N,  and  for  —  =  —  , 

AM  &  C 

PS  ^3 

—  =  3,  and  —  =  100. 

p,  n. 


Q 

SECOND  MODE,  N  =  2 

THIRD  MODE 

N  =  3 

EXACT 

Present 

Method 

Eq.(2.29) 

Kohn,  Krum¬ 
hansl,  and 
Lee,  Ref.  [6] 

EXACT 

Present 

Method 

Eq.(2.29) 

1.0 

1.36 

1.36 

3.  79 

2.  50 

2.  50 

5.  47 

3.  0 

1.  24 

1.  24 

3.  35 

2.  57 

2.  56 

6.  08 

5.0 

1.34 

1.35 

3.41 

2.  51 

2.  49 

6.  22 

FOURTH  MODE,  N  =  4 

FIFTH  MODE,  N  =  5 

Q 

EXACT 

Present 

Method 

Eq.(2.29) 

Kohn,  Krum¬ 
hansl,  and 
Lee,  Ref.  [6] 

EXACT 

Present 

Method 

Eq.(2.29) 

Kohn,  Krum¬ 
hansl,  and 
Lee, Ref.  [6] 

1.0 

3.  77 

3.  76 

6,76 

4.  94 

4.  94 

6.69 

3.  0 

3.  70 

3.  69 

6.  37 

5.  02 

5.  00 

7.  02 

5.  0 

3.  76 

3.  75 

6.  28 

4.  96 

4.  91 

7.  13 

34 


ai 


Q 

EXACT 

Present  Method 

E*  (2.29) 

Kohn,  Krumhansl,  and 
Lee.  Ref.  [6] 

0.27 

0.27 

0.31 

1.91 

1.90 

3.95 

1  0 
,  I.U 

3.48 

3.46 

6.27 

1 

5.26 

5.87 

12.11 

i 

i 

1 

6.58 

7.12 

_ 

13.90 

0.57 

3.68 

6.63 

11.13 


6.78 


7.73 


14.48 


TABLE  HI 


Frequency  Parameter  V  =  atotp/T^)1^2  for  Indicated  Values  of 

b  1  ^3 

Q  =  qa ,  a ,  ri  /n  ,  and  for  —  =  —  ,  —  =  3  ,  and  N  =  5. 
a  i  a  2  p 


FIRST  MODE 

V  =  4 

_ -a  . .  1 _ _ _ 

\!\  =  50 

EXACT 

0.  789 

1.493 

1.  886 

0.  271 

0.497 

0.  609 

\q 

a 

1.  0 

2.  0 

3.  0 

1.  0 

2.  0 

3.  0 

-10.  0 

0.  877 

1.  586 

1.  919 

0.  274 

0.  503 

0.  618 

-  5.0 

0.  879 

1.  590 

1.  923 

0.  274 

0.  504 

0.  619 

0.  0 

0.  929 

1.  650 

1.  983 

0.  292 

0.  530 

0.  639 

0.  5 

0.  936 

1.  719 

2.  101 

0.  279 

0.  532 

0.  677 

1.  0 

1.  811 

3.  622 

5.  4?2 

0.  561 

1.  121 

i.  682 

2.  0 

0.  868 

1.  568 

1.  895 

0  ,73 

0.499 

0.  610 

0.  873 

1.  579 

1.  911 

0.  2  74 

0.  502 

0.  615 
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